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A rigorous three-dimensional relativistic equation satisfied by two-gluon bound states is derived 
from the QCD with massive gluons. With the gluon fields and the quark fields being expanded 
■ in terms of the gluon multipole fields and the spherical Dirac spinors respectively, the equation is 

' well established in the angular momentum representation and hence is much convenient for solving 

the problem of two-gluon glueball spectra. In particular, the interaction kernel in the equation is 
exactly derived and given a closed expression which includes all the interactions taking place in the 
^ I two-gluon glueballs. The kernel contains only a few types of Green's functions and commutators. 

^ ■ Therefore, it is not only easily calculated by the perturbation method, but also provides a suitable 

basis for nonperturbative investigations, 
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I. INTRODUCTION 



It is a prominent feature of Quantum Chromodynamics (QCD) that besides the quark-gluon interaction, there are 
^ : also the interactions between gluons. The self-interaction of gluons suggests that the gluons may form glueballs (the 
bound states of gluons) through their interactions [1,2]. This is an essential prediction which is of decisive significance 
for testing QCD and, therefore, has been arising great interest in searching for the glueballs in experiment [3-7]. But, 
there have been no faithful evidence to be established so far for their existence [8-15]. On the other hand, the property, 
' the mass spectrum, the production and decay of the glueballs have extensively been investigated theoretically. Many 
^p^- approaches were proposed to serve such investigations such as the potential model [16-18], the bag model [19-20], the 
O i' sum rule [21], the Bethe-Salpeter (B-S) equation [22-24] and the lattice simulation [25-30]. However, the theoretical 
D . results given by different approaches are different and even contradictory with each other [31-32]. This situation is 
' attributed mainly to the fact that the quark-gluon confinement has not clearly been understood so far. Commonly, 
^ , it is believed that the lattice gauge approach would give a more reliable prediction because the approach is grounded 
on the first principle of QCD and essentially nonperterbative although certain approximations are inevitably made 
in practical calculations. In addition, it is widely recognized that the B-S equation which is set up on the basis of 
quantum field theory is a rigorous formalism for the bound state problem [33] and suitable to study the glueballs 
" " ' [23-24]. Nevertheless, there are two difficulties in the previous application of this equation. One difficulty arises 
from the interaction kernel in the equation. This kernel was not given a closed expression in the past although the 
expression can be derived by the procedure as demonstrated in a recent publication of one author of this paper [34] . 
Ordinarily, the kernel is defined by a sum of all two-particle irreducible Feynman graphs and can only be calculated 
by the perturbation method. Another difficulty is ascribed to the four-dimensional nature of the equation in which 
the relative time (or the relative energy) is unphysical and would lead to unphysical solutions [35]. So, many efforts 
in the past are paid to recast the equation into a three-dimensional one in the instantaneous approximation[36] or the 
quasipotential approach [37-38]. 

In this paper, we are devoted to deriving a rigorous three-dimensional relativistic equation satisfied by two-gluon 
bound states so as to provide a firm basis for further study. A practical application of this equation to calculate the 
glueball spectrum will be presented in the next paper. The distinctive features of the equation derived are as follows. 
(1) The equation is exactly relativistic, containing all the retardation effect in it, unlike the B-S equation given in the 
instantaneous approximation in which the retardation effect is completely neglected. The equation is derived in the 
equal-time formalism by the consideration that a bound state is space-like and can exists in the equal-time Lorentz 
frame. In this frame, the relativistic equation naturally becomes a three-dimensional one without loss of any rigorism. 
Moreover, different from the B-S equation, the three-dimensional equation derived in this paper is a standard eigen- 
equation of Schrodinger-type. In the position space, it appears to be first-order differential equations. In particular. 
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the interaction kernel in the equation is given a closed expression which is derived by the procedure proposed first in 
Ref. [39] for a two-fermion system and subsequently demonstrated in Ref.[40] for quark- antiquark bound states. The 
kernel derived contains exactly all the interactions taking place in the bound states and is represented in terms of 
only a few types of Green's functions. Such a kernel can not only be easily calculated by the perturbation method, 
but also is suitable for nonperturbative investigations. (2) The three-dimensional equation is established in the 
angular momentum representation. Since a glueball is the state of definite spin and parity, obviously, to investigate 
the glueballs, it is much more convenient to work in the angular momentum representation. In order to express 
the relativistic equation in such a representation, it is necessary to express the QCD in the same representation. 
This can be done by expanding the quantized gluon and quark fields in terms of the gluon multipole fields and the 
spherical Dirac spinors, respectively [41-42]. With these expansions, the vertices in the interaction Hamiltonian can 
all be given explicit and analytical expressions, as will be shown in the next paper. This achievement is due to 
that the integrals containing three and four spherical Bessel functions in the vertices are all calculated analytically 
and expressed explicitly. We would like to note that in comparison with the momentum representation in which 
every glueball state must be separately constructed according to a certain requirement for the Lorentz and OTP 
transformation properties [24] and the vertices would involve the gluon polarization vectors which are not convenient 
to deal with, the angular momentum representation has an advantage that the glueball state can easily be written out 
in a consistent manner and the QCD vertices exhibit the spin structures much clearly. (3) The relativistic equation 
is set up by starting from the QCD with massive gluons. According to the conventional concept of QCD, in order to 
keep the Lagrangian to be gauge- invariant, the gluons must be massless. On the contrary, in most of the previous 
investigations of the glueballs, an effective gluon mass was phcnomcnologically introduced so as to get reasonable 
theoretical results [17,18,23,24]. The gluon mass was supposed to be generated dynamically from the interaction with 
the physical vacuum of the Yang-Mills theory [23] or through strong gluon- binding force [17]. Apparently, these 
arguments would not be considered to be stringent and logically consistent with the concept of the ordinary QCD. 
One of the authors of this paper in his recent article [43] gave a different reasoning that the QCD with massive gluons 
can, actually, be set up on the principle of gauge-invariance without the need of introducing the Higgs mechanism or 
the Stiickelberg fields. The essential points to achieve this conclusion are: (a) The gluon fields must be viewed as a 
constrained system in the whole space of vector potentials and the Lorentz condition, as a necessary constraint, must 
be introduced from the beginning and imposed on the Lagrangian; (b) The gauge-invariance of a gauge field should 
be generally examined from the action of the field other than from the Lagrangian because the action is of more 
fundamental dynamical meaning than the Lagrangian. Particularly, for a constrained system such as the gluon field, 
the gauge-invariance should be seen from its action given in the physical space defined by the Lorentz condition. This 
concept is well-known in Mechanics; (c) In the physical space, only infinitesimal gauge transformations are possibly 
allowed and necessary to be considered. This fact was clarified originally in Ref. [44]. Based on these points of view, it 
is easy to prove that the QCD with massive gluons is gauge-invariant. Moreover, the renormalizability and unitarity 
of the theory have been proved to be no problems [45]. 

The remainder of this paper is arranged as follows. In Section II, the massive QCD and its Lagrangian are 
briefiy described. Section III is used to formulate the angular momentum representation and give the expansions 
for vector, spinor and ghost fields in this representation. In Section IV, the expression of QCD Hamiltonian in the 
angular momentum space will be described and discussed. Section V serves to derive the three-dimensional relativistic 
equation satisfied by the glueballs. In Section VI, we are devoted to derive a closed expression of the interaction kernel 
included in the relativistic equation. In the last section, some remarks will be made. In Appendix, we present a breif 
derivation of the spherical Dirac spinors given in the angular momentum representation. 



II. QCD LAGRANGIAN WITH MASSIVE GLUONS 



In the previous attempt of building up the massive non-Abelian gauge field theory, the following massive Yang- Mills 
Lagrangian density was chosen to be the starting point [43,46]: 

£ = _V^-i.a^ + i^2^aM^a^ (2.1) 

where is the vector potential for a gluon field, 

F^. = - d.Al + gr'"'AlAl (2.2) 

is the field strength tensor in which g is the QCD coupling constant and /"^'^ are the structure constants of color 
SU(3) group and /i is the gluon mass. The first term in the Lagrangian is the ordinary Yang-Mills Lagrangian which 
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is gauge-invariant under the whole Lie group and used to determine the form of interactions among the gluon fields 
themselves. The second term in the Lagrangian is the mass term which is not gauge-invariant and only affects the 
kinematic property of the fields. The above Lagrangian itself was ever considered to give a complete description of 
the massive gauge field dynamics. This consideration is not correct because the Lagrangian is not only not gauge- 
invariant, but also contains redundant unphysical degrees of freedom which must be eliminated by introducing a 
suitable constraint condition. As we know, a massive gauge field has three polarization states which need only three 
spatial components of the four- dimensional vector potential to describe them. In Lorentz-covariant formulation, 
a full vector potential A^^ can be split into two Lorentz-covariant parts: the transverse vector potential Aj, and the 
longitudinal vector potential A^ 

A'' = At^ + Al. (2.3) 

Since the Lorentz-covariant transverse vector potential A"^ contains three- independent spatial components, it is 
sufficient to represent the polarization states of a massive gluon. Whereas, the Lorentz-covariant longitudinal vector 
potential A^' appears to be a redundant unphysical variable which must be constrained by introducing the Lorentz 
condition 



d^Al = 0, (2.4) 

whose solution is 

Al^ = 0. (2.5) 

With this solution, the massive Yang-Mills Lagrangian may be expressed in terms of the independent dynamical 
variables A'^^ 

C = —F^'^'^F^^, + lf.'A^/A^^^, (2.6) 

which gives a complete description of the massive gluon field dynamics. If we want to represent the dynamics in the 
whole space of the full vector potential as described by the massive Yang-Mills Lagrangian in Eq.(2.1), the massive 
gluon field must be treated as a constrained system. In this case, the Lorentz condition in Eq.(2.4), as a constraint, 
is necessarily introduced from the onset and imposed on the Lagrangian in Eq.(2.1) so as to guarantee the redundant 
degrees of freedom to be eliminated from the Lagrangian. Since the action is of more dynamical significance than the 
Lagrangian, the gauge- invariance of QCD should generally be seen from the action given by the Lagrangian in Eq.(2.6) 
or the Lagrangian in Eq.(2.1) constrained by the Lorentz condition in Eq.(2.4). Under the gauge transformations 

sa; = Dfe\ (2.7) 

where 

Df = 6^^% - gr'"'A'^. (2.8) 

Noticing the identity fabcj^anj^b^ _ jg g^gy prove that the action given by the Lagrangian in Eq.(2.1) and 

constrained by he Lorentz condition in Eq.(2.4) is gauge- invariant, 



5S = j d'^xSjC = j d^xO'^d^'AI = 0. (2.9) 

This suggests that the QCD with massive gluons may also be set up on the basis of gauge-invariance principle. 

Now, let us briefly describe quantization of the QCD with massive gluons. This quantization was carried out by 
different approaches in Ref.[43]. A simpler quantization is performed in the Lagrangian path- integral formalism by 
means of the Lagrange undetermined multiplier method which was shown to be equivalent to the Faddeev-Popov 
approach of quantization [44]. For this quantization, it is convenient to generalize the QCD Lagrangian and the 
Lorentz condition to the following forms: 

Cx = - *.gTM«) - m}V - \f^''''F;, 

+ 1^2^aM^a _ l^(^a)2 (2.10) 
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and 



a^A" + aA" = 0, (2.11) 

where for completeness, the quark fields have been included in the Lagrangian in which 'ip and stand for the quark 
fields, r° are the color matrices and m is the quark mass, A°(a;) are the extra functions which will be identified with 
the Lagrange multipliers and a is an arbitrary constant playing the role of gauge parameter. According to the general 
procedure for constrained systems, the constraint in Eq.(2.11) may be incorporated into the Lagrangian in Eq.(2.10) 
by the Lagrange multiplier method, giving a generalized Lagrangian such that 

-ia(A'^)2+A"(a^y4;i+aA") 
= i;{ir{d^ - igT'^Al) - m}V - \f'^'"'F^, + ^i^^A^'^A; (2.12) 
+A"a''A^ + ia(A")2. 

This Lagrangian is obviously not gauge-invariant. However, for building up a correct gauge field theory, it is necessary 
to require the dynamics of the system, i.e. the action given by the Lagrangian (2.12) to be invariant under the gauge 
transformations denoted in Eqs.(2.7) and (2.8). By this requirement, noticing the identity fabc^anj^b^ _ q ^^^^ 
applying the constraint condition in Eq.(2.11), we find 

5Sx = / d^xd''A'^{x)d^'{V''J'{x)9\x)) = 0, (2.13) 

where 

V';f{x)=S-'^d-; + D';f{x), (2.14) 

in which = a/j.^ and D'^{x) was defined in Eq.(2.8). From equation (2.11) we see ^d" A'^ = —X^ ^ 0. Therefore, 
to ensure the action to be gauge-invariant, the following constraint condition on the gauge group is necessary to be 
required 

d^,{Vf{x)e\x))=Q. (2.15) 

These are the coupled equations satisfied by the parametric functions 6°'{x) of the gauge group. Since the Jacobian 
is not singular 

detM ^ 0, (2.16) 

where 

5{d^^Vl^{x)6'^{x)) 

^ ^^'^^ = W{y) l^=° 

= 5''\U^ + a^)5\x -y)- gr''d^{A;{x)S'{x - y)), (2.17) 

the above equations are solvable and would give a set of solutions which express the parametric functions 6'" (x) as 
functionals of the vector potentials A'^{x). The constraint conditions in equation (2.15) may also be incorporated into 
the Lagrangian (2.12) by the Lagrange undetermined multiplier method. In doing this, it is convenient, as usually 
done, to introduce ghost field variables C"'{x) in such a fashion 

e^ix) = ^Cix), (2.18) 

where (, is an infinitesimal Grassmann's number. In accordance with Eq.(2.18), the constraint condition in Eq.(2.15) 
can be rewritten as 



d'^iV^'C) = 0, (2.19) 
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where the number ^ has been dropped. This constraint condition usuaUy is called ghost equation. When the condition 
in Eq.(2.19) is incorporated into the Lagrangian in Eq.(2.12) by the Lagrange multiplier method, we obtain a more 
generalized Lagrangian as follows 

+X''d^Al + ^aiX'^f + Cd^iVfC''), (2.20) 

where C°'{x), acting as Lagrange undetermined multipliers, are the new scalar variables conjugate to the ghost 

variables C°'{x). 

At present, we are ready to formulate the quantization of the QCD with massive gluons. As we learn from the 
Lagrange undetermined multiplier method, the dynamical and constrained variables as well as the Lagrange multipliers 
in the Lagrangian (2.20) can all be treated as free ones, varying arbitrarily. Therefore, we are allowed to use this kind 
of Lagrangian to construct the generating functional of Green's functions 

= ^J D{Al,4,,i,,C\C'',y)exp{i j d^x[Cx{x)+r^{x)Al{x) 

+ V^ry + 7?V + f (a;)C"(x) + C%x)e{x)\}, (2.21) 

where _D(^°, • • • , A°) denotes the functional integration measure, J'^^rj, rj, ^ and are the external sources coupled 
to the gluon, quark and ghost fields and iV is a normalization constant. Looking at the expression of the Lagrangian 
(2.20), we see, the integral over A" (a;) is of Gaussian-type. Upon completing the calculation of this integral, we finally 
arrive at [43] 

= D{Al,ij,i;,C'^,C%)exp{i J d^x[Ceff{x) 

+r'^{x)A;ix) +i^v + vi^ + T{x)c'^{x) + c'^{x)C{x)]}, (2.22) 

where 

£e// = ^{ii^id^, - igT'^A;) - m}v - \f'^'"'f;, + ^i^^a'^'^a; 

-l.^d'^Alf- di^CvfC^ (2.23) 

is the effective Lagrangian given in the general gauges. 

Prom the generating functional shown in Eqs.(2.22) and ( 2.23), one may derive the free gluon propagator as follows 

It is emphasized that when the gluon mass tends to zero, this propagator together with the effective Lagrangian in 
Eq.(2.23) and the generating functional in Eq.(2.22) all immediately go over to the results given in the QCD with 
massless gluons. 

We would like to point out that the above propagator may also be derived from the Lagrangian in Eq.(2.23) by the 
method of canonical quantization [46]. In doing this, we need to use the Fourier representation of the free gluon field 
operator [46] 

A^(-) = j (^{^ ^ k )e-^'^ +e';{k )ar ( k )e^'=1 

^ -[a^o{t)e-'~^'' + a^+{t)e^^''], (2.25) 



2wo{ k ) M 



where k = {ko, k ) and k = (fco, k ) with ko = = \J k + jj? and fco = i^o( k ) = + a^, e^(li) are the 

polarization vectors satisfying the transversity condition 
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fc''e^(fc) = 0, (2.26) 



which means that only three spatial components oi€^{k) are independent and the creation and annihilation operators 
a^( A; ) and a^"^( k) (A = 0, 1, 2, 3) are subject to the following commutation relations 

[a',(l?),a^,(fc^)] = [af (t),a^+(fc^)] = 0. (2.27) 
By making use of the expression in Eq.(2.25) and the above commutation relations, it is not difficult to derive 

iD'^.iix -y) = {0 |T{A^(x)A^(y)}| O) = 1 ^ii^^. (fc)e-''=(--^), (2.28) 

where iD'^{k) is just as that written in Eq.(2.24). 

The propagator in Eq.(2.24) is written in arbitrary gauges. It has been proved that the S-matrices given by the QCD 

with massive gluons is independent of the gauge parameter a [45]. For example, in the tree diagram approximation, 
noticing the transversity condition denoted in Eq.(2.26) and the on-shell property of the gluon states, it is easy to 
verify that the a-dependent term proportional to k^k^, in Eq.(2.24) gives no contribution to the S-matrix elements. 
This fact was early pointed out in Ref.[17]. In view of this fact, we may simply take the Feynman gauge {a = 1) in 
practical calculations. In this gauge, the effective Lagrangian in Eq.(2.23) becomes 

£ = £o + JO.i, (2.29) 

where 

Co = ^Alin + M^)A"^ + ^{i-f'^d^ - m)i, + C"(n + m^)C° (2.30) 

and 

Ci = -i (a^A^ - d^A^A^i'A'''' - ^g'^f°-^'=f'"^^A^''A'"'A'lAl 

+-iij'y^'T"il>A^^ + gf''^''d''C°'C''Al (2.31) 

are the free and interaction parts of the Lagrangian respectively. Correspondingly, the gluon propagator in Eq.(2.24) 
is reduced to 



III. EXPRESSIONS OF THE QCD FIELDS IN ANGULAR MOMENTUM REPRESENTATION 

The fields in the Lagrangian in Eqs.(2.30) and (2.31) may be expressed in the angular momentum representation 

in terms of the eigenfunctions of total angular momenta for the vector, spinor and scalar fields. These eigenfunctions 
are described below. For the vector field, the complete set of the eigenfunctions were already found in the litera- 

ture[41,42]. They include the scalar multipole field A^j^{k'3^) and the vectorial multipole fields ~Aj^{k1^), ^j^(fcl?') 

and ~Ajj^{klt). They are displayed in the following. 

A^jMik-^) = \ -k3j{kr)YjM{x), (3.1) 



Ajj^j{k'f) — —j====LA^jj^{k'^) 



\l^kjj{kry^ jjM 



(£), (3.2) 
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/'2 k . 
- ^2 J _^ ^ ['/jjj+i{kr)Yjj+iM{x) (3.3) 

jj-im{x)], 

Ajj^{k'^) = —VAjj^{k'^) 

= -iJ- jL^ [VJ+Tjj+i{kr)tjj+iM{x) (3.4) 
V IT V2 J + 1 

+\/jjj-i{kr)^ jj-im{x)], 

where J, M mark the total angular momentum and its third component of the vector field eigenfunction, It is the 

momentum of a single particle, k ^ \k\, ji{kr) is the Z-th spherical Bcsscl function with r = L = —iHf x V is 
the orbital angular momentum operator and ^ jimix) denotes the vectorial spherical harmonic function with total 
angular momentum JM and orbital angular momentum Z = J— 1,J, J+1. This function is defined as 

^jim{x) = Cf^yUx)^,, (3.5) 

mq 

where Cf^^ are the Clebsch-Gordan (C-G) coefficients, Yimi^) with x = (6*, </>) are the eigenfunction of orbital 

angular momentum operator L and is the eigenfunction of the spin operator for a vector particle . Since the 
function defined by 

tjM{t)=3j{kr)tjiM{x) (3.6) 

and 

A'>jM{-t)=jj{kr)YjM{x) (3.7) 
respectively satisfy the following equations of motion derived from the first term of the Lagrangian in Eq.(2.30) 

(□ + = (3.8) 

and 

(□ + /z2)A0 = (3.9) 

as we see, the functions in Eqs.(3.1)-(3.4) completely describe the eigenstates of the total angular momentum for a 
free massive gluon field (or for a massless gluon field when the mass /U = 0). Furthermore, as the functions Ajj,^{klt) 
and Ajj^{klt) satisfy the transversity condition 

V-^ = (3.10) 

and are related to each other by 

V X I'^Mik'^) = kljMik^), (3.11) 

the ~Ajn^{k1if) is usually called transverse magnetic multipole field and ~Ajj^{k'x') transverse electric multipole field. 
While the field ~Aj]^(klt) is referred to as longitudinal multipole field because it obeys the condition of longitudinal 
fields V X j4 = 0. It is easy to verify that these vectorial eigenfunctions meet the orthogonality relation 

d-'x'l^jM{k^)^j'M'ik^) = S^x'^ik- k')Sjj'5j^M'^ (3.12) 
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where A, A' = M,E,L label the different modes of the multipole fields. Similarly, for the scalar multipole field, we 

have 

J d'xA^jM{k^)A^j,j^,{k^) = 5{k - k)5jf5^^,. (3.13) 

Therefore, the multipole fields defined in Eqs.(3.1)-(3.4) may suitably be chosen as the basis functions to establish 
the angular momentum representation for the gluon fields. A gluon field operator A^(a;) may be expanded as 

^'(^) = E /'^i[a}l,(A;):I^M(fcl^)e-^-* + a}^M(A;):?jM(fclS^)e^'^*], (3.14) 
where \ = M,E, L, 

Ki^) =12 r ^[a}V(A;)ASM(fc^)e-^"* +ayM(fc)^}M(fc^)e^'^*]. (3.15) 



JM 

The creation and annihilation operators in the above expansions satisfy the following commutation relations 

cA n,\ „c'A'+ 



[a}l^(fc),a},V'(fc')] = 5{k - k')5ec'5xx'5jj'SMM', 

[aS^M(fc),a}?M'(fc')] = [a}1^(fc),a}>M'(fc')] = 0, (3.16) 



[affM{k),a.'yfM,{k')] = S{k - k')Scc'Sjj'6MM", 

[aSV(fc), a}>MKfc')] = [^7M, a}>M. {k')] = 0. (3.17) 
From the transformation of a gluon field under the space inversion [46] 

PX\t, -t)P-'^ = -X\t, -^),PAl{t, lt)P-'^ = Alit, -It), (3.18) 
where P is the space inversion operator and the parity of the multipole fields 

A'jMi-klt) = i-iyA^Mik^), '^jMi-klt) = {-iy+-^Tj^{klt), (3.19) 

where tta = if A = M and tta = 1 if A = i, one may find from Eqs.(3.14) and (3.15) the parities of the annihilation 

operators 

Par/M(fc)P-^ - (-l)V/M(fc).Pa}lf(fc)P-^ = (-l)^+i+-^a}^M(fc) (3.20) 

and the same parity for the creation operators. In addition, from the charge conjugation of the gluon field 
CA^(i, x^)C^^ = —A'jj^{t, ~3f'), we obtain a minus C-parity for the operators aj\j{k) and a'f'j^{k) 

(^^'jM{k)C ^ = —a!jy{k),Ca!f'j^{k)C ^ = — aj'^(fc), (3-21) 

where A = 5, M, L. 

Now let us turn to the angular momentum representation of spinor fields. As shown in Appendix, this representation 
may be set up by means of the spherical Dirac spinors. These spinors are shown in the following. 



v^jMiP^) = (-1)^+-+^ I V '^riPr)^J-Mr-, I ^ ^3_23) 



where JM denote the total angular momentum and its third component of a free fermion, (7 = ±1, e = ^/p"^ + m? is 
the energy of the fermion in which m is the mass and ~]t is the momentum, p = iljj^{x) is the spherical harmonic 
spinor defined as 
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J +a(M-l)+l 



and Uj{pr) is defined by 



^jm(^)= \ — : r , •'-2.--2^ ^ (3.24) 



<(pr) = i-^-t W -pjj_M. (3.25) 



The spinors u'^jj^(p'3^) and Vj;^{plt) are tlie eigenstates of totai angular momentum for a fermion. They respectively 
obey the free Dirac equations of positive energy state and negative energy state. Based on the above expressions of 
the spherical Dirac spinors, it is easy to prove the following orthonormality and completeness relations 



/ d^xu'}'lj{plt)u'},j^^,{p'lt) = S{p - p')S J j'Saa'SMM' 



(3.26) 



J d^xv'}lj{plj^)v'}',^j,{p'li^) = 5{p-p')5jy5„„,5MM', 
J d^xufM{p-xt)v'}',M,{p'^) = J d^xv'}ljip-t)u'}',j^,{p'Tf ) = 0, 
E /o°° dp[u-jj^{p-t)u-jUp-^') + v-jj^{p-^)v-jUp-^')] = 5^{^ - ^')- 

JMcr 

Clearly, the free quark field operators may be expanded in terms of the spherical Dirac spinors 

/•OO 

^{x)=Y. dp[h7M{p)u''jM{p-^)e-'''+d7M{p)v''jM{p-^)e'% (3.27) 

saJM 

^(x)=^/ dp[h7M{p)^JM{p-^)e'''+d7M{pWjM{p-^)e-'% (3.28) 

where s = (c, /), c and / are the color and flavor indices for a quark,It5jy^(p'^) = u'j;^{p'3^)^^Q,v"ji^^{plt) = 
Uj;i^(p"^)+7o and bj]i^ (p), (p) and bj^(p),dj^(p) are the creation and annihilation operators. These oper- 
ators satisfy the following anticommutation relations 



{bSV(p)>byMt(p')} = 5{p-p')5^^,5,,.5jj,5mm', 
{dSV(p),d}',^Mt(p')} = 5{p-p')6^^.5,,,5jj,5mm', 



(3.29) 



with the other anticommutators being zero. 

From the space inversion transformation of the quark field [46] 

P^(^, t)P-^ = r]p-f°il){-lt, t) (3.30) 

and the relations 

7^5m(-P^) = (-1)-'-5«5m(p^),7°^7m(-P^) 

= {-iy-i+\'^j^{p-t), (3.31) 

it is easy to find from Eqs.(3.27) and (3.28) the parity of the quark operators such that 



p^t{p)p-^ = vp{-lY--^^^;M, 

Pdli(p)P-^ = r;p(-l)'^+f +Mf^(p), 

pdt^\p)p-' = M-^y+^+'df^^p), 



(3.32) 



where the phase factor r/p usually is chosen to be r/p = 1. By noticing the charge conjugation transformations [46] 

Bij;{x)B-^ = vcCi^ix), (3.33) 

where C = 17^7*^ and 

v^jm(p'^) = Cu'}Mip-^f,UjMiP^) = Cv-jMiP^f, (3-34) 
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the C-parities of the creation and annihilation operators are easily found from Eqs.(3.27) and (3.28) 

Bbt^\p)B-^ = ^cd.t^^ (p).Bdt^\p)B-^ = vcbt\p). 



(3.35) 



For the ghost fields, considering their scalar character, we can write their expressions in the angular momentum 
representation as follows 



°° dk 

[c'}M{k)A}M{k^)e-'^' + d'}M{k)A^jM{k^)e'"'], (3.36) 



C"{x) = / -=[d'ij^{k)A^j^{k^)e-'-' + c5M(fc)^7M(fc^)e''^*]- (3-37) 
The anticommutation relations for the ghost particle operators in the above are 

{c}jv^(fc),c}'4,(A-')} = 6aa'SjrSMM'S{k - k'), 

{<rjM{k),di+j,{k')} = daa'SjJ'6MM'6{k - k'), (3.38) 

with the other anticommutators being vanishing. 

IV. QCD HAMILTONIAN IN ANGULAR MOMENTUM REPRESENTATION 

In this section, we arc devoted to discussing the QCD Hamiltonian in the angular momentum representation. By 
virtue of the expansions given in Eqs.(3.14), (3.15), (3.27), (3.28), (3.36) and (3.37), it is not difficult to formulate 
the Hamiltonian derived from the Lorentz-covariant Lagrangian written in Eqs.(2.30) and (2.31). However, due to 
presence of the field A^{x), the expression of the Hamiltonian is rather complicated. Considering that the field A^{x) 
is unphysical and will eventually be eliminated from the S-matrix and the B-S equation by the ghost fields as ensured 
by the unitarity of the theory [45] , for the sake of simplifying the representation of Hamiltonian, we would rather to 
start with the Lagrangian obtained from the Lagrangian in Eqs.(2.30) and (2.31) by setting the unphysical fields 

[x) , C {x) and C" {x) to vanish, as was commonly done in the lattice gauge calculations [25-30] and similarly done in 
Ref.[23] where a massive QCD Lagrangian taken in the temporal gauge was used for calculating the glueball spectrum 
in the framework of B-S equation. The Lagrangian we start with is 

C = Co+Ci, (4.1) 

where 

£o = —Ati^ + ii^)A1 + ijii^d^ - m)V 

= ^{Atf + ^^^(V^ - n')At + i^ii^d^ - m)i; (4.2) 



and 

1 

T 

-\9^r'"r'''A'iA]AfA] - ^^iT'^Mf. (4.3) 



jOi = -T^ar'^diA] - djA1)A\A'^^ 



By making use of the canonical variables conjugate to A^, ^ and \E' which are defined by 

Iil = -^ = Al, 
dAl 



= ^ = 0, (4.4) 
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we can write the Hamiltonian density as 



n = iiiAi ^ -L 

= Wo + Uu (4.5) 

where 

Wo = \{Alf - \a1{V^ - ii^)A- + • V + /3m)^ (4.6) 

and 

Hi = X I^) • (V X I^) 

+ -g^f^befcde(ji X A^) ■ (P X A^) (4.7) 

8 

+9^"^ ■ A^ip. 

The first two terms in Eq.(4.6) denote the energy of free gluon fields. The second term gives the energy of free quark 
fields. In the interaction Hamiltonian density (4.7), the first, second and third terms represent the gluon three-line 
vertex, the gluon four-line vertex and quark-gluon vertex respectively. It should be noted that the above Hamiltonian 
includes not only the transverse field A t . but also the longitudinal field A ^ whose presence can be seen from the 
constraint condition in Eq.(2.11) because in the case of A'^ = and a = 1, the condition becomes V-^^ = —A ^ O.This 
is different from the theory proposed in Ref. [23] where only the transverse field exists. 

Now let us derive the expression of the Hamiltonian given by the Hamiltonian density shown above in the 
angular momentum representation. First, we derive the expression of the free Hamiltonian. On substituting 
Eqs.(3.14), (3.15), (3.27), (3.28) into Eq.(4.6), one may get 

Ho = j d^xHo 

= dkoj(t)a.fM{k)a.fM{k) (4.8) 

cXJM 

/>oo 

+ ^ / dpeitWjM {p)bfM{p) + dfM (p)drM(p)]- 

salm 

In later derivations, it is convenient to introduce compact notations for various indices so as to simplify the expres- 
sions of the operators and formulas. For the gluon field, we define 

„cA«.,^ r a5,V(fc) if^ = l 

a.M(fc)-| a}^^(fc)if^ = -l ' ^"-^^ 



XUk^) = \ iiMik^]''^ = ' ^^^^^ 

[ AjMikt) ife = -l 
then, the commutators in Eqs.(3.16) and (3.17) can be unified to write as 

[aj1S(fc)>4M''(fc')] = 5{k - k')dec'6xx'6jj,6MM' sin[(e - O^]. (4.11) 

Furthermore, we define a = {c, X, J, M,k,£^) in which the ^ will be written as later on. In this notation, the 
commutator in Eq.(4.11) will be represented as 

[aa,a^] = A„/3, (4.12) 

where 

Aa/3= d{k-k')5cc'Sxx'Sjj'SMM'sm[{^-Oj] (4.13) 
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and the expansion in Eq.(3.14) becomes 

A^{x) = J2^" 1a{klt)e-'^''"-\ (4.14) 



where 



POO 

^ = ^ / dk, uia = \Jk'^ 

a cXlm •'^ 



(4.15) 



For the quark field, with the definition 



X(fcx^) = -^I^Mik-^)- (4.16) 



= { ;/ J" ' . (4.17) 



drr(p) i/r?=-i 

'^')/ ^^ / "rm(P"^) '7= 1 ^4 iQ\ 
the anticommutators in Eq.(3.29) becomes 

{c7J(p),c;;';j+(p')} = ^(p-p')5..'5..'^jj" W' V- (4.i9) 

When we set a = {s,a,l,m,p,r]) in which the rj will be written as r]a later on, the above anticommutator can be 
simply written as 

{c„,c+} = <5c./3, (4.20) 

where 

Sa/3 = S{p- p)5gg'Saa'Sjj,'6MM'Snn', (4.21) 

and the expansions in Eqs.(3.27) and (3.28) will be represented by 

V'(a;) = I]caWc<(p^)e"'^"^"*, 

a 

where 

J2=J2 dp, eo, = Vp2 + m2. (4.23) 

In view of the notations defined above, the free Hamiltonian will compactly be rewritten as follows 

= H° + ir°, (4.24) 

^s = ^E^"/5:a„a;3:, (4.25) 

a/3 

:4c„ :, (4.26) 

a 

where the symbol :: denotes the normal product and 

= uj{k)5{k-k')5cc'h\'5jj'5MM'{^-5ii')- (4.27) 
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Next, wc derive the expression of the interaction Haniiltonian in the angular momentum representation. In doing 
this, It is helpful to use the transformations for an operator ( a Hamiltonian, a field operator, a creation or annihilation 
operator) among the Heisenberg picture, the Schrodinger picture and the interaction picture 

F„{t) = e'"*Fse-'"* = e'"*e-"'>*Fi{t)e'"'>*e-'"\ (4.28) 

where FH(t),Fs and Fj{t) stand for the operators given in the the Heisenberg, Schrodinger and interaction pictures. 
From Eq.(4.28), it is clearly seen that Fg = Fh{0) = F/(0) which is independent of time. According to the relation 
in Eq.(4.28), we only need to write out the interaction Hamiltonian in the Schrodinger picture (or in the interaction 
picture). In this picture, on inserting the expressions in Eqs.(4.14) and (4.22) into Eq.(4.7), we obtain 



Hi{0) = J d^xni{lt,0) 

= ^ A{af3'y) : 3.^3.03^ : 



J2 B{aPj5) : a^a^sa^a^ : (4.29) 

+ ^C(a/37) : c+c^a^ :, 



where 



A{aPj) = -^-r"" j d?x{X X I^) • (V X X), 



(4.30) 



B{ap-i5) = ^g^r'^r'"' j d^'xiX x I^) • (X; x I^) (4.31) 

and 

C(a/37) = g j d^xw^T'^^Wf) ■ H,. (4.32) 

By making use of the expressions of Aa^Wa and wp, Completing the integrals over it and applying the formulas of 
angular momentum couplings, one may derive explicit expressions of the coefficients A{a[3^), B{a[3^5) and C{aP^). 
The expressions of A{al3"f) and B{aj3^5) will be given in the next paper. As for the coefficient C(a/37),its expression 
will be presented later in other publication. 



V. THREE-DIMENSIONAL RELATIVISTIC EQUATION FOR TWO GLUON BOUND STATES 

The aim of this section is to derive a rigorous three-dimensional relativistic equation satisfied by two-gluon bound 
states. This equation may be derived from the equation satisfied by the following gluon two-time Green's functions 
defined in the Heisenberg picture. These Green's functions include 

G(a+/3+;7-r;ti -t2) = (0+|T{a„(ti)a;3(fi)a+(t2)a+(t2)}|0-), 
G(a-/3-;7+5+;ti -i2) = (0+|T{a+(ti)a+(ii)a^(t2)a5(i2)}iO-), 

G(a-/3+; 7-5+; ii - ^2) = (0+|T{a+(ti)a;3(ti)a+(t2)a5(i2)}|0-), (5.1) 
G(Q+/3-;7+5-;ti -t2) = (0+|r{a,(ti)a+(ti)a^(t2)a+(t2)}|0-), 

where T symbolizes the time-ordered product and |0^) denote the physical vacuum states defined in Heisenberg 
picture. It should be noted that the creation and annihilation operators in Eq.(5.1) are those defined in Eq.(3.15). 
For later convenience, we will use the notation defined in Eq.(4.9) for the gluon operators. With this notation, instead 
of the above Green's functions, we may start, in a consistent way, from the following Green's function defined by 

C?(a/3;7^;ti -i2) = (0+|T{a„(ii)a^(ti)a^(i2)a5(i2)}|0-), (5.2) 
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where 



where be a creation operator or an annihilation one. 

On differentiating Eq.(5.2) with respect to time ti^ it is found that 

J J_G(a/3; 7,5; h - t2) = iSih - t2)S{a(3, jS) 

+(0+m[(^^a„(ti))a^(ii) (5.3) 
+aa(ii)(iJ:a;3(ii))]a^(i2)a5(t2)}|0-), 

S{ap,'y6) = {O+\[aa{h)a0{h),a^{t2)aLs{t2)]\O-) 

= (O+|[aaa0,a^a5]|O-) ^^■^> 

is time-independent owing to the time displacement shown in Eq.(4.28) and the restriction of 5{ti — t^)- According 
to the relation in Eq.(4.28), we have 

i^a„(t) = e'^*[aa,/f]e-^^*. (5.5) 

Since in the commutator is independent of time, we may use the Hamiltonian given in Eqs.(4.24)-(4.26) and (4.29) 
and the commutation relation in Eq.(4.12) to compute the commutator [aa,i?]. When the result is substituted in 
Eq.(5.5), we obtain 

i-^^aiti) =Y,'^pa^aa3ip{tl) + E h{P<^'r)^ctr ■ ap(tl)a^(fi) : 
pa par 

+ E f2{p(TTX)Aa\ : ap{ti)a„{ti)ar{t-i) : g^ 

parX ^ ^ 

ptTT 

where Aaa was defined in Eq.(4.13) and the coefficients fi are defined by 

/i(pCTr) = A{paT) + A{pt(t) + A{Tpa), 

/2(parA) = BiparX) + B{pa\T) + B{pXaT) + B{XpaT), (5.7) 
/3(pa-T) = C{paT). 

Clearly, the expression of i^a^(fi) may be written out from Eq.(5.6) by replacing a with (3. Inserting this expression 
and that given in Eq.(5.6) into Eq.(5.3), one may find an equation of motion obeyed by the Green's function denoted 
in Eq.(5.2) such that 

i J_G(a/3; 7,5; h - h) - i6(h - t2)S{a(3, j6) 

+ E l^paAai3,a\G{pX; jS; ti - t2) 

pcrX 

+ E .9i("/3;pcrA)Gi(pcrA;7(5;ti - ia) ,^ 

paX 

+ E 92{aP\ P<yTX)G2{p<TTX;-i5;ti - t2) 

parX 

+ E 53 (a/?; paX)G3{paX; j5;ti-t2), 

paX 

where 

Aa/3,aA = AaaSffX + A/S^SaX, (5.9) 

gi{aP;paX) = E /i(P^'r) Aa/3,rA, 

r 

g2{aP] parX) = E f2{p'JT0)Aa0,ex, (5.IO) 
gi{a(j; paX) = E /3(po-'r)Aa/3,rA, 

T 

and the Green functions Gj are defined by 

Gi{paX\'i5;ti - 12) 

= (0+|T{: a^(ii)a<,(ii) : ax{h)a^{t2)a5{t2)}\Q-), (5.11) 
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G2{p(TT0;^5\tx - t2) 

= (0+|T{: a^(ti)a^(ti)a^(ti) : ae(ti)a^(i2)a5(t2)}|0-), 



(5.12) 



(0+|T{: c+(ii)c<,(ti) : aA(ti)a^(i2)a5(t2)}|0-). (5.13) 



By the Fourier transformation 

f +00 



1 /■"^°° 

Gi(a, 5;io) = - dte'^'Gi{a, 5; t), (5.14) 

* J —00 

the equation in Eq.(5.8) will be of the form in the energy representation 

+ J2 (^pa^ap.a\G{pX;-fS;uj) 

+ E 5i("/3;pcrA)Gi(pcrA;7(5;tj) 

pctA (5.15) 
+ J2 g2iaP; pcrTX)G2ip<7TX;^S;u!) 

parX 

+ J2 gi{al3] paX)G?,{paX;^6]uj). 

Let us look at the second term on the right hand side (RHS) of the above equation. From the definitions of Wp^ and 
Aaff^crX given in Eqs.(4.27), (5.9) and (4.13), we see, only when ^ to pa and A^^^ are nonvanishing. When we 
define 



U!oc = { 

the second term mentioned above gives 



X] ^pcrAai3,axG{pX; ^5;ui) = (wc + L0f3)G{aP; 7(5; w). (5.17) 

paX 

Thus, Eq.(5.15) can be written as 

(u) -u)a- W/3)G'(a/3; jS; iv) 

= S{a/3,j6) + J2 gi{al3; paX)Gi{paX;j6;u) 

paX 

+ E 92iaf3: paTX)G2{p<JTX;-fS;uj) (5.18) 

parX 

+ E £/3(a/3;P<^A)G3(pcTA;7(5;a;). 

paX 

In order to obtain a closed equation, it is necessary (actually possible) to introduce effective interaction kernels Aj 
in such a fashion 



^ gi {aP; paX)Gi {paX; 7^; w) 



pcrX 



^Ai(a/3;p(7;w)G(/9c7;7(5;w), (5.19) 



g2{af3; paTX)G2{p(TTX;'yS;uj) 

parX 

Y,Mal3;pa;uj)G{pa;jd;uj), (5.20) 



15 



pcrX 

= AsM; per; to)G{pa; jS; u). (5.21) 

pa 

With the kernels defined above, Eq.(5.18) becomes 

(w - - uji3)G{al3; 7(5; w) 
= SiaP, 7<5) + P^; ^)Gipa; ^S; to), (5.22) 

pa 

this is just the equation satisfied by the Green's function G{ap; 'y6; co) in which 

3 

K{al3; pa; w) = ^ Ai{ap; pa; w). (5.23) 
j=i 

is the total interaction kernel. 

The equation obeyed by glueball states may be derived from the above equation with the aid of the Lehmann rep- 
resentation of the Green's function. Suppose {|n)} form a complete set of glueball states, noticing the transformation 
in Eq.(4.28) and H\n) = i?„|n), we can write 

G{aP;^d;h-t2) 

= EiHh - t2)(0+|a«a^|n)(n|a^a5|0-)e-^^''(*^-*=) (5 24) 

+0{t2 - ti)(0+Ka5|n)(n|a„a^|0-)e-*^''(*=-*i)}. 

Substituting in Eq.(5.24) the representation of the step function 

and then performing a Fourier transformation with respect to time, we can get from Eq.(5.24) the Lehmann repre- 
sentation such that 

G{aP;,S;u;) = ^[^1!^^^ _ X.s{n)x.,{n) 

where En is the total energy of a glueball state and 

Xa/3(") = (0+|aaa;3|n) , x„^(n) = (n|aaa^|0-) (5.27) 

are the B-S amplitudes describing the two-gluon bound states. On replacing the G{aj3;^5;uj) in Eq.(5.22) with its 
Lehmann representation, then multiplying the both sides of Eq.(5.22) with u — En and finally taking the limit 
CO — > En, we eventually arrive at 

{En -iOoi- (^0)Xa0{n) = Y P^''' En)Xpa{n). (5.28) 

pa 

This is just the wanted three-dimensional relativistic equation for the two-gluon bound states. It should be pointed 
out that the above equation actually represents a set of coupled equations which may be written out by setting ^q, in 
the index a and in the index [5 to be itl.Each of the equations is manifestly of Schrodingcr-typc, i.e., a standard 
eigen-equation. In the next paper, these coupled equations will be reduced to an equivalent equation satisfied only 
by the B-S amplitude for which the two gluons are in the positive energy states. 

VI. CLOSED EXPRESSION OF THE INTERACTION KERNEL 

In this section, we are devoted to derive the effective interaction kernel appearing in Eq.(5.28) and defined in 
Eq.(5.23), giving a closed expression of it. For this derivation, we need equations of motion which describe the 
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where 



evolution of the Green's function G{a[3\^5]ti — ^2) and those shown in Eqs.(5.11)-(5.13) with time t2- Taking the 
derivative of the Green's function in Eq.(5.2) with respect to we have 

i^^GiaP; 7(5; ti - ts) = -iS{h - t2)S{aP, 7^) . . 

+(O+|T{a0(ii)aMii)[(*J7a^(i2))a^fe)+a^fe)*J;a,(t2)]}|O-), > 

pa 

+ "£ /2(p(ttA)A^a : a,(i2)a„(t2)a.(i2) : ^^'^^ 

parX 

+ E Mp'^rJA^r ■■ cj{t2Kit2) :, 

par 

here the coefficients are the same as represented in Eqs.(5.7). The expression of i-^3is{t2) can directly be written 
out from Eq.(6.2) by replacing 7 by 5. When this expression and the one given above are inserted into Eq.(6.1), we 
are led to 

t^G{a(};^6;ti - fs) = -iS{h - t2)Siap,jS) 
+ E ^pa^tS,a\G{a0; pX; h - t2) 

per A 

+ J2 Gi{aP;paX;ti - t2)gi{pa\;'i5) 

paX (6-3) 
+ G2{a(3;paTX;ti - t2)g2{p(TTX;'y6) 

parX 

+ E Gs{af3;paX;ti - t2)gs{p(jX;^5), 

paX 

where S'(a/3, 7*5) was defined in Eq.(5.4), the coefficients Qi are the same as those represented in Eq.(5.10) since the 
following equality holds 

A; 7^) = 5i(7(5; p, A) (6.4) 
and the Green functions Gj are defined as follows 

Gi{a(3; paX;ti - 12) 

= (0+|T{a„(ii)a^(ti) : a^(i2)a<,(t2) : aA(i2)}|0-), 

G2{aP; parX; ti - 12) ,„ 
= (0+|T{a„(ti)a;3(ti) : a^(i2)a,(t2)a,(t2) : aA(i2)}|0-), ^^-^^ 

Gs,{al3] paX]ti - ^2) 

= (0+|T{a„(ii)a/3(ti) : c+(i2)c,(t2) : aA(i2)}|0-). 
By the Fourier transformation denoted in Eq.(5.14) and noticing Eq.(5.17), we obtain 

{ui -{-ui^ +u}s)G{al3;^5;Lo) = S{aP,jS) 

- E Gi{al3; p(7X;uj)gi{p(7X;jS) 

paX 

- G2{af3; p(jrX:uj)g2{po-rX;j5) (6.6) 

parX 

- E GaiaP; paX;uj)g3{p(7X;'-fd). 

paX 

According to the same procedure as formulated in Eqs.(6.1)-(6.6), one can derive the equations of motion obeyed by 
the Green's functions denoted in Eqs.(5.11)-(5.13). Because the same expressions of the differentials i-^a^{t2) and 
i-^as{t2) are employed in all the derivations, the equations of motion for those Green's functions may immediately 
be written down by referencing the equation in Eq.(6.6). The equation for the Green's function in Eq.(5.11) is 

{uj + uj^ + u)s)Gi {paX; jS; u)) = Si [paX, jS) 

- E Gii{paX; p,VT;uj)gi{p.UT;^5) 

- E GuipcrX; pi^TK;uj)g2{pJ^TK;j6) (6-7) 

- E Gis{p(jX\ pvT\uj)gz{p,VT;-i5), 

fJ-VT 
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(6.9) 



where 

5i(paA,7<5) = (0+|[: a^a,, : aA, a^aaJlO"), (6.8) 

Gii{paX; fii'T;u!),Gi2{p(TX; ij,i'tk;uj) and Gi3{paX; fiVTiw) are the Fourier transforms of the following Green's func- 
tions 

Gii{paX;fj,vT;ti - ^2) 

= (0+|T{: a^(ti)a^(ti) : Six{ti) : a^(t2)a.(t2) : a^(i2)}|0-), 
Gi2(pcrA; fXUTK; ti — ^2) 

= (0+|T{: RpitMh) : axih) : a^(i2)a.(i2)a^(i2) : a«(i2)}|0-), 
Gi3{paX;piyT;ti - 12) 

= (0+|T{: ap(ti)a<,(ti) : axih) : c+(t2)c,(i2) : a^(t2)}|0-). 
For the the Green's function in Eq.(5.12), we can write 

{(jj + co^ + ijjs)G2{p(7tX; ^6; ui) = S2{p<jtX, j6) 

- Yl G2i{p(ttX; piy9;u))gi{piy9;'y6) 
live 

- G22{pcrTX; iivK6;uj)g2{lJiVK6;^5) (6.10) 

- J2 G2z{p(JtX; pv9;uj)gz{pv0;-i5), 

where 

S2{paTX,'yS) = (0+1 [: apa^^a,- : aA,a^a5]|0") (6.11) 
and the Green's functions on the RHS of Eq.(6.10) are defined as 
G21 (pcrA; pi^6; ti — ^2) 

= (0+|T{: ap(ti)a^(ii)a,(ti) : aA(ti) : a^(i2)a,(t2) : ae(t2)}|0-), 

G22(/5o-tA; iivkB] ti - ^2) 

= (0+|r{: ap(ti)a,(fi)a,(ti) : aA(ti) : a^(f2)a,(t2)aK(t2) : ae(f2)}|0->, 
G23.{p<ttX] p,v0]ti - 12) 

= {0+\T{: ap{ti)a„{h)ar{h) : aA(ti) : c+{t2)c,{t2) : a0{t2)}\O-). 



(6.12) 



For the Green's function in Eq.(5.13), we have the equation 

{uj + Uj + us)G3{paX;"fS;uj) = SsipaXj^S) 

- J2 G3i{paX; iJLVT;oj)gi{nvT;^5) 

- G32{p<TX;pi'TK;uj)g2{niyTK;j6) (6.13) 

- E G33{p(jX; p.vT;uj)g3{p.vT;^5), 

where 

S3{paX,^5) = (0+1 [: c+c, : aA,a^a5]|0-) (6.14) 
and the Green's functions on the RHS of Eq.(6.13) are defined by 

G^iipaX; puT-^ti - 1,2) 

= (0+|r{: c+(ti)c,(ti) : axih) : ap(t2)a.(t2) : a^(t2)}|0-), 

G32ip<TX;piyTK;ti-t2) . . 

= (0+lr{: c+{tMh) : aA(ti) : a^(i2)a.(t2)a,(t2) : a«(t2)}|0-), ^^"^^^ 
G33{paX;puT;ti - ^2) 

= (0+|T{: c+{ti)c^{h) : aA(ti) : c+(t2)c,{t2) : a^(t2)}|0-). 

Now we are ready to derive the interaction kernels. Multiplying the both sides of Eq.(5.19) with + cos) and 
then applying Eqs.(6.6) and (6.7), we have 
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Ai (a/3; pa; u) [S {pa, 7(5) - ^ Gi {pa; p,v\; w)gi {p,vX; j6) 

per fivX 

= Y9AotP'-: P^>'^)[Si{P'^>^^l^) -YGii{paX; iJ.PT;u))gi{ni^T;^S) (6.16) 

paX HVT 

- ^ Gi2{paX;iJ,i'TK;(j)g2{lJ.i'TK;jS) 

- X] *^13 (po-A; /iz/r; a;)s'3 {p,VT; 7^)] . 

In order to obtain the expression of the kernel Ai (a/3; pa; w) from the above equation, it is necessary to eUminate 
the kernel in the second, third and fourth terms on the left hand side (LHS) of the above equation. Considering that 
the Green's function G{ap;'y6;w),as a matrix G, has an inverse G~^{a(3;'y6;u)), then from Eq.(5.19) we can get 

Ai(a/3;7^;a;) = ^^5i(a/3;p<TA)Gi(paA;/iz/;w)G"^(/iz/;7(5;a;). (6.17) 

pcrX 

Substituting Eq.(6.17) into the second, third and fourth terms on the LHS of Eq.(6.16), then moving these terms to 
the RHS of Eq.(6.16) and finally acting on the both sides of that equation with the inverse of the matrix S{pa,^S) 
which is assumed to exist, we eventually arrive at 

Ai (a/3; jd; u) = A^^^ (a/3; jS; u) + A^^ (a/3; 7^ ; u) 

+A^^\a/3;-fd;u), (6.18) 



where 



and 



AW(a/3;7^;a;) = ^^5i(a/3;paA)5i(paA,/.z/)5-i(Mz/,7'^), (6-19) 

pcrX p-v 

Af)(a/3;7<5;a;) 
= - X] 51 (a/3; PC^A){^ Gil (pcrA; iiut; uj)gi {jiur; On) 

paX 6v fifT 

+ G 12 {pa \; pi'TK;(j)g2{p-i'TK;9Tr) (6.20) 

+ GisipcrX; pur; Lu)gz{iiiyT] Oi^)}S~^{Ott, 7(5), 



Kf\al3:-i8;u) 

^ ^ ^ ^ 31 (a/3; pcrA)Gi {paX; pv; oj)G~^{pv; et; u) 

pcrX p.f ei 6v 

^{J2^^ '^^91 iCvC; On) 

+ J2 G2{ei;^vCx;oj)g2{^vCx;0n) (6.21) 
+ J2 G3{ei; ^vC; ^)93{^lC, 07r)}S-\97T, jS). 
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By the same method as described above, the kernel A2(q!/3; u)) can be derived from the relation in Eq.(5.20) by 
using the equations in Eq.(6.6) and (6.10). The result is as follows 



where 



A2(Q/3;7<^;a;) = A«(a/3;7^;a;)+A(')(a/3;75;a;)+A(3)(Q/3;75;a;), (6.22) 
A«(a/J;7(5;a;) = ^ ^52(a/?;p<7TA)52(parA;/xi/;a;)5-i(Mi/,7^), (6.23) 

parX fif 

= - ^ ff2(a/3; po-TA){y~^ G21 (po-rA; ni^e; oj)gi {live; 0tt) 

pa-rX 9tv pve 

+ G2i{p<yT\; jjLVSx; uj)g2{lJ.ys>c; On) (6-24) 
-^G23{p(JT\;iiv£;Lo)gz{liv£;e'K)}S~^(6'K,'^5), 



and 



pcrrX pv ei Sir 

^ iYl '^1 '^^91 i^vC, Ott) 

+ Y.G2 {et; inCx, ^)52 (€^Cx; Ot^) (6.25) 

+ E ^3(ei; er?C; u^)ff3(?r?C; ^7r)}5-i(0^, 7^). 
(.vC 

Analogously, the kernel h.^{a(i;^5;w) may be derived from Eqs.(5.21), (6.6) and (6.13). The result is shown below 

A3 (a/3; 7(5; w) 

= A«(a/?;75;a;) + Af (a/?;75;a;) + Af (a/3;75;a;), (6.26) 



where 



aW (a/3; 7(5; a;) = ^^53(a/3;paA)53(paA,H5-i(/zi.,7(5), (6.27) 

paX pv 

Af (a/3;7^;a;) 
= - X] ^93{ali\p(T\)(^G3i{pa\\pvT\uj)gi{p.vT;e-K) 

prrrX Stt A'f'' 

+ ^ G32{pcrX:, pvTK;u)g2{piyTK;9Tr) (6.28) 
+ X] ^33 (po-A; /iz/r; a;)s'3 (/iz/r; 6'7r) } 5"^ (^Itt, 7^) , 

Af (a/3;7,5;c.O 
= XI XI XI XI ^3 (a/?; P<jX)G3 {paX; jiv; uj)G~^{nv] et; u) 

pcrX pv ei 6v 
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^ iYl ^^91 {^VC, Ott) 

+ ^ G2 {ec; ^vCx; ^)92 {^vCx; On) (6.29) 
+ ^ G3(et; ^r?C; a;)53(^??C; ^7r)}5-i (^TT, 7^). 

Based on the expressions described in Eqs.(6.18)-(6.29), the total kernel can easily be written out. Using the matrix 
notation for the Green's functions and the other functions, the kernel is represented as 

3 



33 3 

= {E - J2 9iGii9j + Yl giGiG-^Gjgj}S-\ (6.30) 

i=l iij=l i)J=l 

This is just the closed expression of the interac;tion kernel. According to the general argument as presented in 
Refs. [33,34,39,40], the third term in the above expression plays the role of eliminating all the B-S reducible (two- 
particle reducible) diagrams contained in the first two terms. In fact, the relations in Eqs.(5.19)-(5.21) and the 
following ones 

Gjgj = GAj (6.31) 
are inserted into the last term in Eq.(6.30), it is seen that 

3 3 

QiGiG-'Gjgj = Y ^iG^j = KGK, (6.32) 

which exhibits the typical structure of the B-S reducible part of the interaction kernel. Therefore, the kernel shown in 
Eq.(6.30) is truly B-S irreducible, consistent with the conventional concept. The equation in Eq.(5.28) and the kernel 
in Eq.(6.30) will be employed, in the next paper, to calculate the glueball spectrum in the ladder approximation. 



VII. CONCLUDING REMARKS 



In this paper, the exact three-dimensional relativistic equation for two gluon glueball states and its interaction 
kernel have been derived from the QCD with massive gluons. When the gluon mass tends to zero, the equation and 
its kernel will naturally go over to the ones for the QCD with massless gluons. As shown in Eq.(5.28), The equation 
derived is a standard eigenvalue equation of Schrodinger-type. In the position space it appears to be a set of first-order 
differential equations. This kind of equation was given for fcrmion systems [39,40], but never formulated for boson 
systems in the past. It should be noted that since the equation is derived in a special equal-time Lorentz frame, it is 
certainly not Lorentz-covariant even though the equation is rigorous and includes all the retardation effect in it. In 
order to derive a Lorentz-covariant equation, one may start from the four-time Green's function 

Ga0js{T,t;T',t') = (0+|T{a„(ii)a/3(t2)a^(i3)a5(t4)}|0-), (7.1) 

where 

T' =^{t3+U),t' = t3-U. (7.2) 

Differentiating the above Green's function with respect to ti and t2 and employing the expression of the differential 
i-^aa{ti) as shown in Eq.(5.6) and the similar expression for i^a^(i2), one may obtain two equations for the 
Green's function Ga()'ys{T,t;T' ^t'). Adding the both equations and subtracting one equation from another, following 
the same procedure as formulated in Sec.V, it is easy to derive two equations satisfied by the B-S amplitude Xai3{T, t) 
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which respectively describe the evolutions of the state with respect to the center of mass time T and the relative time 
t. In the equal-time frame, owing to the relative time being absent, we are left only with the equation with respect to 
the center of mass time as given in Sec.V. It is interesting to note that either the equation with respect to T or the 
equation with respect to t, appears to be a first order differential equation of Schrodinger-type in the position space 
whose solutions are determined merely by the initial condition of the B-S amplitudes at the time origin. This is an 
essential feature of the equations mentioned above which is different from the B-S equation. The latter equation is a 
higher order differential equation and hence, like the Klein-Gordon equation, has unphysical solutions with negative 
norm as pointed out in the previous literature [35]. This is because the solutions of B-S equation are determined 
not only by the initial amplitudes at the time origin, but also by the time-differentials of the amplitudes at the time 
origin. 

Another point we would like to note is that unlike the Dyson-Schwinger equation, the relativistic equation derived 
in this paper is of a closed form. In particular, the interaction kernel in the equation is given a closed expression. The 
expression contains only a few types of Green's functions and vacuum expectation values of the operator commutators. 
They are unambiguously defined in the Heisenberg picture and each of them can independently be calculated by the 
perturbation method without concerning other Green's functions. Especially, the kernel represents all the interactions 
taking place in the bound states and, therefore, are suitable for nonperturbative investigations because the Green's 
functions and vacuum expectation values of the commutators, in principle, are able to be evaluated by a certain 
nonperturbative method as suggested by the lattice gauge approach. 

At last, it would be pointed out that although the equation in Sec.V and the kernel in Sec.VI are derived in the 
angular momentum representation, they suit to formulate the equation and the kernel in the momentum representation 
as long as the angular momentum quantum numbers in the indices a, /3, • • • are replaced by the momentum ones. 
That is to say, the equation and the kernel formally remain unchanged in the both representations. 
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IX. APPENDIX: SPHERICAL SPINORS IN THE ANGULAR MOMENTUM REPRESENTATION 

In this appendix, we intend to give a derivation of the spherical Dirac spinors which are used, as basis functions, 
to establish the angular momentum representation for fermion fields. It is well-known that in the relativistic case, 
unlike the helicity, the spin of a free fermion is not a good quantum number. However, the total angular momentum 
operator of the fermion commutes with the Hamiltonian. Therefore, it is meaningful to discuss eigenfunctions of the 
total angular momentum which satisfy Dirac equation. Let us start from the positive energy spinor Ua{'^) which is 
the solution to the Dirac equation {id^p^ — m)u{'J^) = 0. This spinor is taken to be [46] 




which is normalized in such a fashion: u+(~j?)us(~p^) = 1. The negative energy spinor can be given by the charge 
conjugation Ws(l?) = CusCj^)^ . Suppose Yim{p) with p = |"^| = {9, ip) and ips are the orbital angular momentum 
and the spin eigenfunctions respectively; the total angular momentum eigenfunctions may be constructed in the 
momentum space by the C-G coupling 

^'jM iP) = E ^lJ:isYlm {p)Vs , (A.2) 
Is 

where s = §,(t = ±1,Z = J± 2. Noticing the representation 

n=(;)'^-i=(;) (A-3) 

and employing the explicit expressions of the C-G coupling coefficients for different values of a which may be found 
in the textbook, one may derive from (A.2) the expression as follows 



22 



J+a(M-l) + l 



Y 2J-<7+l ^J-f,M+iW 



By making use of the Dirac spinor in (A.l) and the eigenfunctions represented in (A. 2) or in (A. 4) with respect to 
the momentum we may construct the spherical Dirac spinor in the position space through the following Fourier 
transformation 



Substituting the expansion 



e^t-^ =4nJ2 i'ji iPr)Yi*m {x) (A.6) 



Im 



and the expressions written in (A.l) and (A. 4) into (A.5), considering 

^008 6* sin^e" 
, sin 616^^ - cos 61 

and 



^ ( cos 9 sin 9e ^'^ \ / a r,\ 



^■pn'}M{P) = -^jL{P), (A.8) 

which is easily proved by utilizing the familiar recursion formulas for the spherical harmornic functions, it is not 
difficult to derive the expression shown in Eq.(3.22). The expression of the function Vjj^{p'3f) may be derived by 
the charge conjugation denoted in Eq.(3.34). The result was written in Eq.(3.23). It would be noted that the 
eigcnfunction ilj]^j{x) in Eq.(3.24) which is defined in the position space and appears in Eqs.(3.22) and (3.23) is of 
the same form as the function njj^^{p) in (A. 4) which is defined in the momentum space. 
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